Dynamics of topological defects in a spiral: a scenario for the spin-glass phase of 

cuprates 
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We propose that the dissipative dynamics of topological defects in a spiral state is responsible for 
the transport properties in the spin-glass phase of cuprates. Using the collective-coordinate method, 
we show that topological defects are coupled to a bath of magnetic excitations. By integrating out 
the bath degrees of freedom, we find that the dynamical properties of the topological defects are 
dissipative. The calculated damping matrix is related to the in-plane resistivity, which exhibits an 
anisotropy and linear temperature dependence in agreement with experimental data. 
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The discovery of four static incommensurate (IC) 
charge and spin peaks by neutron scattering experiments 
in Nd-doped La2- K Sr K Cu04 (LSCO) confirmed the 
proposal that in these materials the holes form verti- 
cal/horizontal charge stripes |2j, breaking the paradigm 
of homogeneous charge distribution. Moreover, the ob- 
servation of IC dynamical spin correlations in supercon- 
ducting LSCO (x > 0.05) at the same wave vectors as ob- 
served in the Nd-doped compounds 0] placed the stripe 
picture on a more solid footing and raised further ques- 
tions concerning the coexistence of magnetism and super- 
conductivity. Later, two IC elastic magnetic peaks have 
been observed within the spin-glass (SG) phase of LSCO 
(0.02 < x < 0.05) H, and their interpretation in terms of 
diagonal stripe formation appeared to be rather natural, 
given that the value of the magnetic IC follows the same 
linear dependence on doping as observed within the su- 
perconducting regime Q. However, no charge order has 
ever been measured in the SG regime, and the picture of 
well-ordered stripes, very far apart from each other due 
to the very low doping concentration, seems quite im- 
probable because disorder and frustration effects, which 
act to destabilize the stripes, are expected to dominate 
in this phase. 

A different explanation of the two IC peaks in the SG 
phase is that these may arise from the formation of a spi- 
ral magnetic order which breaks the translational sym- 
metry in the spin but not in the charge sector |5j. In 
this model, the randomly distributed holes act as frustra- 
tion centers for the underlying anti-ferromagnetic (AF) 
background, generating a dipole moment. A fraction of 
these dipoles may order ferromagnetically, while the oth- 
ers may remain disordered. If the number of ordered 
dipoles increases linearly with the doping concentration, 
which would appear plausible, this model can describe 
the linear variation of the IC magnetic wave vector with 
doping 4] , without invoking charge ordering [5j . 

Recently, transport properties along both diagonal- 
directions in the Cu lattice have been investigated in 
high-quality, detwinned LSCO samples, and a rather 



strong anisotropy was detected jfjjf . However, this experi- 
ment provides no definitive evidence in favor of the stripe 
picture, because, as we will show, a similar anisotropic 
behavior is expected for transport parallel or perpendic- 
ular to the spiral axis. 

In this Letter we investigate the formation of topolog- 
ical defects in the spiral state that couple to excitations 
of the magnetic environment and diffuse across the sys- 
tem. We assume that the charge carriers are attached to 
these defects and therefore the description of defect dy- 
namics corresponds ultimately to the electrical transport 
properties of cuprates in the SG phase. We evaluate the 
damping matrix, which is related to the mobility of the 
topological defect, and show that its behavior at high 
temperatures follows the experimentally observed tem- 
perature (T) dependence of the in-plane resistivity. In 
addition, we account for the anisotropy in the spin-wave 
velocity in the directions parallel and perpendicular to 
the spiral axis, and find that, in agreement with experi- 
ments, it leads to an anisotropy in the resistivity. 

The magnetic properties of charge transfer insulators, 
such as undoped La2Cu04, are usually studied within 
the quantum non- linear sigma model (NLtrM), which 
can correctly represent the long wavelength modes of the 
Heisenberg model. This continuum model describes slow 
fluctuations of the locally well defined staggered mag- 
netization n (with n 2 — 1). When holes are added to 
the system, the magnetism becomes more complex be- 
cause those act simultaneously as a source of disorder 
and dipolarlike frustration. The dipolar frustration can 
be described via a coupling of the dipoles to the gradient 
of the order parameter n of the NLerM. By accounting for 
the random distribution of dipolar centers, the resulting 
Hamiltonian favors the formation of a spiral phase, with 
a non-zero average twist c^n of the AF order parameter 
and a concomitant alignment of some of the dipoles [f| . 

The spiral ground-state breaks the 0(3) spin sym- 
metry completely. Its order parameter, which belongs 
to the SO (3) group, is given in terms of an orthonor- 



mal basis nfc, with k — 1,2,3 and n\n a q 



= 5. 



kq 



2 



The latter can be related to an element g of SU(2) 
through Ufc = (l/2)tr {a a ga k g^ x ), where a a are Pauli 
matrices. It is also convenient to introduce the fields 
A c ^ = (l/2i)tr (a a g~ 1 d u g) , which are related to the first 



derivatives of through d^n^ 



Here 



/i stands for one-time and two-spatial coordinates. The 
Lagrangian describing the spiral state is given by 



L = Af I dx\\dx 



(as-. 



a 2 _ 2*2 

_A ± C||A|| 



(1) 



where M = 2JS 2 /c\\c±, J is the AF exchange and cj_, en 
are the spin wave velocities perpendicular and parallel to 
the spiral axis, with (c||/cj_) 2 = cos (Qa) ||. Here a is 
the lattice constant and Q is the IC magnetic wave vec- 
tor observed by neutron scattering. Lagrangian may 
be mapped to an isotropic form by introducing the coor- 
dinates r = (x,y) with x — y/cnjc±x±, y 



We then find L = Af J d 2 r (A§ 



v/cl/c||£||. 
c 2 A 2 ), where a = x,y 



and 



is the isotropic spin-wave velocity. This 



C = y/C±C\\ 

rescaling leads to an anisotropic elementary cell, which 
will be neglected in what follows. Further comments on 
the anisotropy effects are included in the conclusions. 

An appropriate treatment of the spiral state has to take 
into account the existence of topological defects. In prin- 
ciple, in (2 + 1) space-time dimensions two types of static 
topological defects can exist: vortices and skyrmions. 
Skyrmionlike excitations are classified according to the 
second homotopy group of the order parameter space, 
which in our case is trivial, iT2{SO(3)) = 0. Thus, there 
are no skyrmions in the spiral state and the only possible 
defects are Z 2 -vortices, given by the first homotopy group 
ir\(SO(3)) = Z2 Q. These defects originate from the chi- 
ral degeneracy of the spiral, which can rotate clock- or 
counter-clock- wise jlOj . Their description is provided by 
the elements g — g s (^)g E (e) of the SU(2) group, with 
g s = cxp (im a a a ^/2) and g £ = exp (ia a e a /2). The pa- 
rameters m}S? are related to the local spin, where m is a 
constant unit vector and if? is a scalar field 0] . Our aim 
is to determine the topological defect solutions *f?2 V of the 
saddle point equations of the spiral, as well as the vec- 
torial field e describing fluctuations around this solution. 
In terms of these fields, the Lagrangian becomes 



L = AT J d 2 r (L + L x ) , L = \(d^f 



(2) 



d »£ + o i 9 ^ x e) 



with dnAd a B = d t Ad t B - c 2 VAVB. This Lagrangian 
can be treated perturbatively because \e\ is small. We 
begin by considering the static limit of Lq, where the 
corresponding equation of motion becomes the Laplace 
equation W 2l f?(r) = 0. A nontrivial solution ^i„(r, R) = 
arctan(j/— Y) / (x—X) has the form of a vortex centered at 
R = (X, Y), and its energy E cx InL/a diverges with the 
system size L. Unbound, free topological defects cannot, 



therefore, exist at low temperatures. However, the two- 
vortices solution ^>2v — ^luO^Ri) — ^it,(r, R 2 ), which 
describes a bound pair of topological defects (vortex and 
antivortex), centered respectively at R x and R 2 , has a 
finite energy E( 1 $> 2 v) oc ln(rf/a), where d = |Ri — R 2 | 
is the modulus of the pair relative coordinate, which we 
take as a constant. Up to an irrelevant constant, the two- 
vortices solution can be expressed in terms of the center 
of mass R = (Ri + R 2 )/2 and relative coordinate d, 
yielding * 2 „(r) = arctan{[(r-R) x d] J[(r-R) 2 -c£/4]}. 

We will apply the collective coordinate method |l2( to 
study the dynamics of a bound pair of topological defects 
in a spiral. The first task is to determine the dynam- 
ics of the fluctuations around the defect. By evaluating 
Li{ 1 $2v) and expressing the vectorial field in terms of po- 
lar coordinates, e = (ecos(9,esin0,£ z ), we find that the 
e z component is free and that 9 = 5 , 2 „/2. The remaining 
equation of motion then reads {df — c 2 [V 2 + l^(r)]} e = 
0, with V(r) = (V* 2l ,) 2 /4. Its solution has the form 
e(r,t) = J2nm Qnm(t) t| nm ( r )i where r] nm are the eigen- 
functions of the operator 



- c 2 [V 2 + V(r)}i] mn = Lul m i] nm . 

The solutions of Eq. have the form rj„ 

Vk nm /SL[ H^\k nm r) + e- 2 ^Hi 2) (k nm r)] 



(3) 



where 



H„ ' 2 ^ are Hankel functions of the first and second kinds, 
d n represents a phase shift for the mode n, and $ is a po- 
lar angle. The excitation with frequency u) nm has the 
dispersion u> nm — c|k nm |. Equation © admits zero- 
frequency modes associated with the continuous trans- 
lational symmetry of the model. A consistent treat- 
ment of them requires the use of collective coordinates. 
Thus, the center of the pair R is promoted to a dy- 
namical variable R(i), namely \&(r) — > ^(r — R(^)) and 
e(r, t) -> e(r - R(t), t) = J2 nm q„ m (i)?7 nm (r -K(t)). 

A lengthy but straightforward procedure [lj| leads to 
the effective Hamiltonian (in units of 2 JS 2 /c 2 ) describing 
the pair of topological defects interacting with the bath 
of magnons, 

nm 

where M — J c? 2 r(V^E' 2l ,) 2 is the mass of the topological 
defect, P and p nm are respectively the momentum con- 
jugate to R and to the coordinate q nm of the excitation 
with frequency u nm , and P E = Y^ nm ,kl Pnm G«m,W Qkl- 
The constants G n m,ki are related to the eigenfunctions 77 
by G nm ,ki — J d 2 rrji c i'Vr] nm . The classical Hamiltonian 
is quantized by introducing the usual creation and an- 
nihilation operators a\ lm and a nm . Restricting the anal- 
ysis to the case of elastic scattering and retaining only 
the terms at quadratic order in , a n 



Am 1 a nm we obtain 



H 



P 2 

2M 



~M 



E 

nrn,kl 



(5) 



3 



nm^nm ^nm 



and D r 



J;l 



where H b = £ nm Hu. 
G 

i i in 

Our aim is to investigate, in the low-energy sector, 
the effective dynamics of the vortex-pair in the presence 
of the excitation bath. We compute the reduced den- 
sity matrix of the pair by integrating out the degrees 
of freedom of the bath of magnons, for which we em- 
ploy the Feynman- Vernon path-integral formalism. The 
evolution of the density matrix for the full system is de- 
scribed by p{t) = exp(— iHt/h)piO)ex.piiHt/h). For the 
sake of simplicity, we use the factorisable initial condi- 
tion p(0) = p v (0) /5b(0), where p v and pb represent, re- 
spectively, the initial vortex-pair and bath density ma- 
trices. This condition implies that the pair of topologi- 
cal defects and the excitations do not interact at t = 0. 
The bath degrees of freedom are supposed to be initially 
in thermal equilibrium, pb(0) — e~P Hb /tt\e~P Hb ], where 
/3 = h/ksT. After evaluating the trace over these, we ob- 
tain the reduced density-matrix operator for the vortex- 
pair 



p(x, y, t) 



J J dx'dy'J(x,y,i;x',y',0)^(x',y',0), 



where the super-propagator J has the form 



J 



/;-/: 

Jx' Jy' 



p ye i[So[x]-So[y]] jP [ X)y ]_ 



(6) 



5o[x] = J di'(M/2)x 2 is the action associated with the 
free motion of the vortices and y] is the influence 
functional, which describes the effect of the excitations 
on the dynamics of the topological defect pair. 

Defining the center of mass and relative coordinates 
v = (x + y)/2, u = x y and using the Born approxima- 
tion, one finds, after a lengthy calculation that 



.F[x, y] = exp exp [$] , 



(7) 



where 



$ = 2V f dt'f dt"e" v (t' -i"K(iV(t"), *= I dt ' I dt"l^it' -t")u^it')u v it"), 

6^(t) = ^ N nmD% mM D^ mtkl $m(LJ nm -UJ k l)t, l> 1V {t) = N nmD^ m kl D^ nhkl COS iuj nm -UJ k i)t. 



nm.kl^O 



The boson occupation number is N nm = (e' 3 "™™ — 
The effective action that describes the motion of the vor- 
tex pair in the presence of excitations may then be ob- 
tained by inserting Eq. Q into the superpropagator 
We find that S e f f — So [x] — So [y] + h& and the corre- 
sponding equations of motion indicate that the dynamics 
of the defects is damped, 

v"(t) + V f dt'Y v {r - if) i v (lf) = 0, 

where the matrix 7 All/ (i) = (2Ti/M)e ,J " l '(t) is a general- 
ization of the damping coefficient. The decaying part 
exp[$] of the influence functional is related to the diffu- 
sive properties of the defect, with the diffusion matrix 
V^it) — fie if). The evaluation of the diffusion is 
analogous to that of the damping matrix, and at low 
T the two parameters are related by the fluctuation- 
dissipation theorem. Here we will concentrate on 7 only. 
We first introduce the scattering function S^iui, u>') = 



EG: 



nm 



kl 



5(u; — u> nm )8{u>' — u>ki), which allows the 



description of the problem in terms of continuous fre- 
quencies. Because the model we study is isotropic, the 
damping matrix is diagonal 7 MI/ (t) = 7(i)(5 M!y , and 7(t) 
is obtained directly from the first derivative e(t). After 
introducing the substitution £ = (cj + lo')/2, ( = lo — u>', 
and considering only processes for which £ << 1, we find 
that 7(t) = j(T)5(t) with 



2nH r At frt^K) 



(9) 



e 



The phase shifts 5 n of the eigenfunctions y n m can be 
calculated using the Fredholm method 0, [Uj. Their 
evaluation requires the computation of the matrix ele- 
ments An = (n\V\n), because S n = arctan[7Tv4 n /(l + 
An)]- Here \n) are the eigenstates of the unperturbed 
operator in Eq. J3J). In the limit of a small vortex- 
antivortex separation, d <C r, the potential reads V{r) — 



4 



d 2 /[4(r 2 + gP/4) 2 ]. One then obtains 

•An(£) = -Trffife [/ n (£d/2c)# n (£d/2 C )] , 

where /„ , K n are the modified Bessel functions of the 
first and the second kinds, respectively. Using the Drude 
model and assuming that the charge carriers are attached 
to the defects, we can relate 7 to the inverse of the mo- 
bility, p^T 1 = Mj/e. We then find 



nh/3 x 
ec 2 ^— ' 

n=l 



(10) 



n+l 



AO 



1+A n + A n+ 1 + (1 + ir 2 )A n An +1 ' 



This expression leads to a vanishing /i" 1 as T — > 
because in this limit there are no magnons to scatter the 
topological defects and the latter behave as free particles. 
However, this limit will never be reached in systems like 
LSCO because at Tj ~ 30K the dipoles freeze and the 
charge becomes localized, leading to an upturn of the 
resistivity. At high T we obtain 



TTk B T 



ecr 



i + BKO'' 



(ii) 



which indicates that the resistivity p = (/in^e) -1 varies 
linearly with T, as observed experimentally 6, 14]. Here 
n n denotes the charge carriers concentration. This be- 
havior is expected to hold for T < Tkt, where Tkt is 
the Kosterlitz-Thouless temperature, at which the vortex 
anti- vortex pairs will eventually unbind. An attempt to 
roughly estimate Tkt has provided Tkt ~ JS 2 0. Us- 
ing that J - O.leV and 5=1/2, we find T KT ~ 300K. 
Despite of the crudeness of the estimates, this result com- 
pares quite well with the experimental data |l4|. which 
shows deviations from linearity at T ~ 400K. Given that 
Tkt represents the temperature at which few pairs start 
to unbind, it is not surprising that the linear behavior 
can hold up to higher T. 

The order of magnitude of the resistivity can also 
be promptly evaluated from our calculations. Although 
the exact microscopic values of the parameters d and c 
are not available, we nevertheless provide an estimate 
of these. Taking hc/ahs ~ Tkt and d ~ a, we find 
pT 1 ~ 0.1Vs/cm 2 for T — 200K, in good agreement with 
Fig. 2 of Ref. 0. 

In order to provide also a comparison with the exper- 
imental data concerning the anisotropy in dc transport, 
we should include the spin-wave anisotropy in the spiral 
state. The analytical calculations which led to Eqs. © 
- <|11|) cannot be performed for an anisotropic magnon 
bath. Nevertheless, we provide a rough estimate of this 
effect by considering the dependence of the inverse mo- 
bility in Eq. I|ll|) on the spin-wave velocity c. Recalling 
that B n is a function of £d/2c, we find that pT x oc cr 1 . 



One could then expect that in the original anisotropic 
system /ij 1 oc cj 1 and pj 1 oc cj 1 . This result can also 
be understood on physical grounds: since dissipation is 
provided by the excitations of the magnon bath, spin- 
waves with higher velocity are less effective in scatter- 
ing the defects. As a consequence, the resistivity in the 
spiral direction is larger than that in the direction per- 
pendicular to it, 7_l/7|| = Px/p\\ — c \\/ c ± — \/cos (Qa). 
For x = 0.04, we find y/cos (Qa) = 0.98, in agreement 
with transport measurements by Ando et al. |6j, which 
found that for 100K< T < 200K the resistivity along 
the a-axis in LSCO is slightly smaller than along the b- 
axis (in orthorhombic coordinates). We note that the IC 
peaks observed in neutron scattering correspond to the 
6-direction, which coincides with the spiral axis because 
the breaking of translational symmetry within the spiral 
picture has its origin in the spiral chirality (which can ro- 
tate clock- or counter-clockwise). Thus, p± = p a , pu = 
and p a < Phi as experimentally observed. This result in- 
dicates that the anisotropy measured in the SG phase 
does not provide evidence for the existence of diagonal 
stripes, but instead, is also the result expected within a 
more realistic, albeit more complex, model which does 
not need to appeal to charge order at such small doping 
concentrations. 

In conclusion, we propose a description of the trans- 
port properties in the SG phase of cuprates based on 
the dissipative dynamics of topological defects in a spiral 
state. Using the collective-coordinate method, we derive 
the effective action of the topological defects, which are 
coupled to a bath of magnons. The scattering of magnons 
by the potential provided by the topological defects leads 
to a dissipative motion for these defects. Under the as- 
sumption that the holes are attached to the defects, the 
corresponding damping matrix is calculated, and is re- 
lated to the in-plane resistivity. Its T dependence and 
anisotropic behavior are in agreement with the available 
experimental data, indicating that further investigations 
are required to distinguish between spiral spin states and 
diagonal stripes in the SG phase of cuprates. 
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